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Two multiscale- type turbulence models are implemented in the PAB3D solver. The models are based on 
modifying the Reynolds-averaged Navier-Stokes equations. The first scheme is a hybrid Reynolds-averaged- 
Navier-Stokes/large-eddy-simulation model using the two-equation ke model with a Reynolds-averaged-Navier- 
Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length 
scale. The second scheme is a modified version of the partially averaged Navier-Stokes model in which the 
unresolved kinetic energy parameter A is allowed to vary as a function of grid spacing and the turbulence length 
scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale 
resolution possible for a given flow on a given grid for partially averaged Navier-Stokes. It has been found that the 
prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f k varies 
between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier-Stokes 
turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, 
and validation examples are presented to demonstrate the enhancement of PAB3D’s time-accurate turbulence 
modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for 
accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier-Stokes 
results and experimental data for high-temperature jet flows. The current results show promise for the capability of 
hybrid Reynolds-averaged Navier-Stokes and large eddy simulation and partially averaged Navier-Stokes in 
simulating such flow phenomena. 


I. Introduction 

T HE limited capability of the Reynolds-averaged Navier-Stokes 
(RANS) approach, combined with eddy-viscosity turbulence 
models to simulate unsteady and complex flows, has been known for 
some time. The RANS assumption is that most of the energy is 
modeled through the turbulence transport equations and is resolved 
in the grid. RANS overpredicts the eddy viscosity, which results in 
excessive damping of unsteady motion. Consequently, the eddy 
viscosity attains unphysically large values due to unresolved scales 
and suppresses most temporal and spatial fluctuations in the resolved 
flowfield. One of the approaches used to overcome this problem is to 
provide a mechanism for the RANS equations to resolve the largest 
scales of motion. Among several methods, the detached eddy 
simulations (DES) [1], the hybrid large eddy simulation (LES) [2,3], 
the limited numerical scheme (LNS) [4], and the partially averaged 
Navier-Stokes (PANS) [5] provide the mechanisms needed to satisfy 
this requirement. 

In an attempt to increase the fidelity and accuracy of the PAB3D 
code^ [6,7], a hybrid turbulence model RANS/LES [2,3] and PANS 
[8] have been added. Abdol-Hamid and Girimaji [8] explored a new 
approach to improve the accuracy and robustness of PANS in 
creating a simulation of an unsteady flowfield. They accomplished 
this through the development and implementation of a two- stage 
procedure to efficiently estimate the level of scale resolution possible 
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for a given flow on a given grid for PANS and other hybrid models. 
This capability was implemented in a general applied aerodynamics 
computational fluid dynamics (CFD) research code, PAB3D, for the 
simulation of unsteady flows. 

Both large temperature and pressure fluctuations have a profound 
effect on turbulence. Although several models have been developed 
to account for the effect of pressure fluctuations (compressibility 
correction models), very little has been done to account for large 
temperature fluctuations. This has led to poor CFD prediction of 
nonisothermal flows. For high-temperature jet flow, the standard 
turbulence models lack the ability to predict the observed increase in 
growth rate of the mixing layer [9,10]. Several researchers [11-16] 
have modified one or more terms of the transport equations to obtain 
better agreement in high- temperature flows. These modifications 
directly or indirectly affect the closure terms of the turbulent heat flux 
pufi and stresses pupTj. Theis and Tam [ 11 ] changed several 
coefficients in the turbulent transport equations. However, such 
extensive modifications of model coefficients completely change the 
characteristics of the equations and may cause deficiencies in flow 
prediction accuracy for other problems. Other attempts to sensitize 
the turbulence model to temperature fluctuations involve more 
sophisticated closure for the turbulent heat flux term appearing in the 
average energy equation [14-16]. Explicit algebraic nonlinear heat 
flux models have also been tested for this purpose. These models 
have been successful in some fully developed high-temperature 
turbulent flows. A simpler approach was to model the value of C jJL as 
a function of the total temperature gradient in the flow [11]. The 
concept behind this approach was to postulate that the large-scale 
density nonuniformity in the flow would introduce local mixing 
instability and added turbulence stresses. The modification in [ 17 ] 
was successful in predicting the hot-jet mixing rate over a wide range 
of temperature ratios between the jet and the ambient air. The authors 


'Data available online at http://www.asm-usa.com/PAB3D [retrieved 
15 October 2007]. 
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of this paper thought that this variable C ^ would not be necessary in 
the multiscale turbulence model approach because the computations 
in the code would automatically take care of the large-scale density 
gradient effect. 


II. Approaches 

The governing equations of the time-averaged formulation 
include the conservation equations for mass, momentum, and energy 
and the equation of state. In the present study, the perfect gas law was 
chosen to represent the air properties, and the eddy -viscosity concept 
was used to model the Reynolds stresses. The mass, momentum, and 
energy conservation equations of the time-averaged equations can be 
written in a conservative form as follows: 
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In the case of RANS equations, a standard turbulence model 
(STM) such as the two-equation ke turbulence model is used: 
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The boundary conditions with n as the wall distance for s and k at 
the wall are 


£ wall — V l 



^wall — 0 


In the present paper, we use the simple eddy diffusivity (SED) 
approach, which is based on the Boussinesq viscosity model. This 
approach is used to model all the scalar diffusion terms appearing in 
the RANS and standard ke equations. For the heat flux term, the SED 
is written as follows: 


— z P y t dT 

pu t e = — 

<r T dx ; 


a T = 0.9 


The turbulent stress components were formulated as 
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For the purpose of this paper, the RANS turbulent viscosity was 
defined as 


..RANS 

v t 


= f c - 


(4) 


A. Two-Stage PANS Approach 

The PANS model [5] was developed to overcome the grid 
dependency associated with the use of other hybrid turbulence 
models (HTM). In its original form, PANS [5] replaced the two- 
equation turbulence model by solving for the unresolved kinetic 
energy k u and the dissipation e u . The k u equation is identical to the 
original k equation. In the s equation (2), the following coefficients 
were used to change the two-equation model to the HTM, which 
becomes known as the PANS formulation through the following 
changes: 

C S 2 — ^jr(C s 2 ~ C e i) + C e i C>k = ~J~ &e = fy (J £ (5) 

J£ J £ J £ 

where 

k p 

r ^_u_ r c pj_ 

J k j J £ 

/c s 

represent the ratios of the unresolved kinetic energy and dissipation 
to the total kinetic energy and dissipation, respectively. It was only 
natural to use f k and f £ to quantify the PANS filter with respect to 
RANS. Therefore, f k and f B are used as the resolution control 
parameters for PANS. The physics of turbulence dictates that 

(DNS) 0 </*</,< 1 (RANS) 


The original formulation [5] used constant values for the 
unresolved kinetic energy parameter f k and unresolved dissipation 
rate parameter f s . The users will select values for these parameters 
and refine the grid until the flow solution converges toward a solution 
target. This could be very time-consuming for resolving complex 
three-dimensional flows. In the present paper, we will discuss an 
approach to define the unresolved kinetic energy parameter. Abdol- 
Hamid and Girimaji [8] introduced a two-stage approach to estimate 
the values of the unresolved kinetic energy parameter. Here, we will 
highlight the basic concepts of this approach. Based on a simple 
dimension analysis, we assume that the turbulent viscosity may be 
related to the total kinetic energy k, s, S, and A as 


which leads to 
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C h is a model coefficient, which needs to be calibrated. In the present 
paper, a value of 1 will be used to evaluate the model. Here, we 
summarize the guidelines to be followed for the sequential two-stage 
procedure. These guidelines are completely dependent upon flow 
complexity. For the first stage, 

1) Complete a three-dimensional or two-dimensional simulation. 

2) Use unsteady or steady calculation; high-order schemes are not 
required. 

3) Choose a desired level of an allowable RANS turbulence model 
(one-equation, two-equation, algebraic stress, full Reynolds stress, 
etc.). 

For the second stage, 

1) Conduct a three-dimensional simulation. 

2) Use unsteady calculation; high-order schemes should be 
considered. 

3) Use hybrid models (DES, Hybrid RANS/LES, PANS, etc.). 

The users need to use the same flow conditions, boundary 

conditions, and grid resolution for both stages of the procedure. 
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B. Hybrid RANS/LES Approach 

Nichols and Nelson gave an example of a hybrid RANS/LES 
turbulence model. This method was implemented in conjunction 
with Menter’s SST two-equation turbulence model and was termed a 
multiscale (MS) model. In the present paper, this hybrid model is 
used with the two-equation model described in Eqs. (2) and (3). The 
turbulent length scale used in this implementation was defined as 


L = max 


/ (vf A™ k V2\ 

( 6 V n ’ —) 




duj du; 


£2 s 

The subgrid turbulent kinetic energy was defined as 

k LES =f d k 

The damping function was defined as 

f d = { 1 + tanh[27r(A - 0.5)]}/2 

where 

A = - 1 


dxj dxj 


' 1 + [/,/A] 4 / 3 1 + X 4/3 

X is the unresolved characteristic ratio, and 
A = max(A x , A y , A z ) 

The eddy viscosity was then calculated from 

v t =f d vf ANS + (l-fM ES 

vp s = min(vf ANS , 0.084 A V k LES ) 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 
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Note that this hybrid model allowed the transition from RANS to 
LES as a function of the local grid spacing and the local turbulent 
length scale predicted by the RANS model rather than as a function of 
the grid spacing alone. This allowed the model to detect whether it 
can resolve the turbulent scales present on the existing grid before its 
transition over to the LES mode. 


f k =l RANS Regions 

10 



10 

f k <l PANS Regions 


a) Unresolved turbulent kinetic energy parameter, 
f k contours 



b) Snapshot of velocity contour 

Fig. 1 Single-nozzle predictions using PANS formulation. 


III. Results and Discussions 

Two test cases were selected to evaluate the present modifications. 
The first test case is a single subsonic jet operated at the design 
pressure condition. We selected this case to compare RANS, PANS, 
and LES approaches in simulating subsonic jet. This case provided 
the calibration of the PANS approach. The second test case is a 
multistream subsonic jet configuration. This case addressed the 
complexity of modeling multistream flows. We used a 0.03 time step 
based on the radius of the nozzle and the freestream acoustic speed. 
We ran 4000 iterations to initiate the unsteadiness of the jet flow. 
Then the solutions were averaged over the last 20,000 iterations. It 
was observed that approximately four subiterations per physical time 
step produced the optimal convergence per iteration. However, the 
physics of the specific problem will dictate the subiteration number 
for other cases. In the present results, four subiterations typically 
reduced the residual by two orders of magnitude at that time level, 
with no improvement using more iterations. 

A. High-Temperature Single Nozzle Jet Flow 

The present study used the benchmark experiments performed by 
Bridges and Brown [18] at the NASA Glenn Laboratory with the 
flow condition as indicated in Table 1 for the core and freestream. 


Table 1 Experimental subsonic condition 



T t , °R 

P„ psi 

M 

Core 

1400 

17.68 

0.55 

Freestream 

540 

14.3 

0.01 


DATA 


PANS 



RANS 


Fig. 2 Two-dimensional time-averaged velocity contour results from 
RANS, and PANS formulations compared with experimental data. 


The computational mesh is a full three-dimensional grid with 120 
cells in the circumferential direction. The computational domain is 
divided into 48 blocks. The superfine mesh had a total of 4,000,000 
cells. Grid points are clustered near the solid surfaces and around the 
shear layer. We used a uniform stream wise grid spacing for 
2 < x/Dj < 12. The value of y + for the first cell off the surface 
varied between 0.2 and 2. During the course of simulating this case, 
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x/Dj 

Fig. 3 Normalized jet centerline velocity results using RANS, RANS/ 
LES, and PANS formulations compared with experimental data. 



X/Dj 

Fig. 5 Normalized unresolved turbulent kinetic energy k u compared 
with total kinetic energy k t using PANS. 



■ 

uu/U 2 Experimental Data 

▲ 

w/Uj Experimental Data 

• 

ww/U 2 Experimental Data 

- uu/U 2 PANS 

- w/Uj PANS 



- ww/U 2 PANS 



Fig. 4 Normalized jet centerline turbulent kinetic energy results using Fig. 6 Normalized normal stress components ( uu , vv , and ww) 
RANS, RANS/LES, and PANS formulations compared with compared experimental data, 
experimental data. 


we tried single- and double-precision calculations. We found no 
significant differences in the results. Also, we tried 4, 8, and 12 
subiterations for the dual-step time-accurate approach. We also 
found no significant differences in the results. Based on these results, 
we will use single precision and four subiterations in all of the 
presented results in this paper. 

First, we used the RANS formulation to get time-averaged 
quantities to calculate the characteristic length-scale ratio. This ratio 
varies in space and is used to produce the unresolved kinetic energy 
parameter Figure la shows the distribution of this function. This 
parameter identifies the RANS and PANS regions. The RANS 
regions are defined with the parameter set at a value of 1. The PANS 
regions are the remaining flow domains in which f k values are less 
than one. We use this parameter in solving the PANS formulation. 
We calibrate the medium grid (1,000,000 cells) to get the velocity 
profile to closely resemble the experimental data for the velocity 
distribution. We found that C h should be 1.05 for the present test 
case. All the calculations performed hereafter used the same value of 
C h = 1.05. 

Figure lb shows the snapshot of a two-dimensional contour on the 
X-Y plane for the velocity u/Uj using the PANS formulation. The 


result shows the unsteady behavior of the jet flow as it interacts with 
the external flow. A similar observation was found as the RANS/LES 
formulation was used. In the RANS prediction, there was not a 
significant difference between the snapshot and the time-averaged 
flow quantities. In this case, the snapshot solution is similar to that 
shown in Fig. 2. This was caused by the fact that RANS overpredicts 
the eddy viscosity, resulting in excessive damping of unsteady 
motion. Consequently, the eddy viscosity attains unphysically large 
values due to unresolved scales, and suppresses most temporal and 
spatial fluctuations in the resolved flowfield. 

Figure 2 shows the two-dimensional slice of the time-averaged 
results of using RANS and PANS formulations compared with 
experimental data. There were no significant differences between the 
time-averaged and the snapshot of the unsteady RANS result. The 
PANS results and similarly RANS/LES (not shown) produce a good 
comparison with experimental data. Figure 3 shows a comparison of 
the jet centerline velocity using RANS, RANS/LES, and PANS 
formulations compared with experimental data. Both RANS/LES 
and PANS produced better results compared with the RANS 
solution. Figure 4 shows the jet centerline turbulent kinetic energy. 
Both PANS and RANS/LES overpredict the total turbulent kinetic 
energy. We define total kinetic energy k t as 
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Table 2 Experimental subsonic condition- 



T t ,° R 

P „ psi 

M 

Core 

1498 

21.72 

0.8 

Fan 

647 

24.36 

0.9 

Freestream 

530 

14.7 

0.28 


a Please compute the Mach numbers for the core and fan. 


k t = k r + k u = \[uiUi — Ui Ui\ + k u 

In the case of the RANS results, the resolved kinetic energy k r is 
zero. Both RANS/LES and PANS provided a mechanism for the 
RANS equations to resolve the largest scales of motion. Figure 5 
shows the reduction of the unresolved kinetic energy k u compared 
with the total value. In general, the RANS produces the same values 
for all three normal stress components. Figure 6 shows the 
comparisons and distributions of the normal stress components ( uu , 
vv, and ww ) using PANS. In this case, PANS produced different 
values for the three normal stress components with a trend similar to 


the experimental data. The levels of the vv and ww components are 
very close to the experimental data. The predicted uu component is 
much higher than the experimental data. DeBonis [19] reported the 
similar observation that the LES simulation produced a much higher 
value of uu compared with M = 0.9 hot-jet experimental data. 

We attempted to provide a grid convergence study for the single- 
jet case. We tried three grid levels: 1,000,000 cells (medium), 
2,000,000 cells (fine), and 4,000,000 cells (superfine). In this case, 
the grid was refined in the flow direction while holding the filter 
coefficient C h [see Eq. (8)] constant. This causes the filter width to be 
refined along with the grid. The smallest value of f k was around 0.1. 
Under such a procedure, the solution converged to a DNS because 
the filter width approaches zero along with the grid spacing. As a 
result, this approach attempts to calculate different turbulent scales 
explicitly, whether or not the numerical resolution was sufficient to 
accurately compute those scales. The result of this study was shown 
in Fig. 6. The medium and fine grids provided reasonable grid 
convergence. However, the superfine grid completely diverged from 
the other solutions, as shown in Fig. 7. This was due to limitations of 



Fig. 7 Normalized jet centerline velocity results using PANS 
formulation compared with experimental data. 


Experimental Data 
PAB3D PANS Fine 
PAB3D PANS Superfine (1) 



°0 5 10 15 20 

X/Dj 

Fig. 8 Normalized jet centerline velocity results using PANS 
formulation compared with experimental data. 



Fig. 9 Configuration 1 baseline round core nozzle with fan nozzle. 
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X/D 


Fig. 11 Comparisons of computed normalized centerline stagnation 
temperature and data [10]. 

the numerical scheme to handle such flow physics. In such a case, a 
higher-order scheme could be more effective. Now we freeze the 
value of f k from the fine-grid solution used for the superfine grid. 
Figure 8 shows the comparisons of this approach with fine and 
superfine grids. The result of this approach improves the superfine- 


grid solution and now it matches the fine-grid solution. This 
equivalent is to use f £ of value 0.85 instead of 1.0. 


B. High-Temperature Multistream Nozzle Jet Flow 

The second test configuration includes a separate fan and core 
nozzle flows at a bypass ratio of five with an external plug. One set of 
data was selected from the reported test results [20], with the flow 
condition as indicated in Table 2 for the core, fan, and freestream. 

This test configuration was part of a comprehensive investigation 
on jet exhaust noise due to the pylon-chevron-jet interaction [10], 
which was tested at NASA. As previously discussed, the present 
PANS approach was tuned on the axisymmetric grid to the round 
nozzle experimental results. For the limited space and scope of this 
paper, we will discuss the comparison with only the axisymmetric 
configuration shown in Fig. 9. 

The computational domain for the solution extended from x/D c = 
—6.3 to 31.6 in the axial direction and 63 D c in the radial direction, 
where D c is the diameter of the baseline core nozzle, 12.80 cm. The 
origin, x/D c = 0.0, was set at the exit of the fan nozzle so that the exit 
of the core nozzle was at about x/D c = 0.5. The computational mesh 
is a three-dimensional- shaped grid with 120 cells in the 
circumferential direction. The computational domain is divided into 
92 blocks. The mesh had a total of 7,750,000 cells. Grid points are 
clustered near the solid surfaces and around the shear layer. The 
value of y + for the first cell off the surface varied between 0.16 and 
1 . 8 . 

Computational solutions were obtained for RANS, hybrid RANS/ 
LES, and PANS formulations. The simulated conditions were set to 



Fig. 12 Comparison of the stagnation temperature prediction with data [10]. 
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coordinate with the data presented in [10]. The computational results 
were computed for a freestream Mach number of 0.28. For the PANS 
calculations C h = 1.05 and C h = 0.95 to study the sensitivity of the 
prediction to the variation of C h . 

Comparisons of computed stagnation temperatures at the 
symmetry plane show that both RANS/LES and PANS produced 
faster mixing than the RANS solution, as shown in Fig. 10. 
Comparisons between computed stagnation temperature and 
experimental data are shown in Figs. LL and 12. The RANS solution 
produces slower mixing than the experimental data. On the other 
hand, RANS/LES produces better agreement with the data than the 
RANS results. In Fig. 11, the calculations were compared with the 
experimental centerline values. We also found that the PANS 
solutions are sensitive to the variation of C h . As shown in Fig. 12 , 
RANS/LES was able to more accurately predict the temperature 
flowfield and produce a result closer to the experimental data, 
whereas the RANS approach overpredicted the temperature in the 
core region. The results from PANS were in reasonably good 
agreement with experimental data. 

IV. Conclusions 

The hybrid RANS/LES and PANS turbulence models are 
relatively new and will need to be exercised for a wide variety of 
problems to determine their accuracy before they become an 
accepted tool for fluid dynamics modelers. They seem to offer much 
for unsteady flow applications, but issues such as grid sensitivity 
need to be further addressed. It is hoped that more effort will go into 
these models in the near future, so that they can be matured for use in 
everyday applications. The new capabilities have the potential to 
improve the accuracy and robustness of creating a simulation of an 
unsteady flowfield. This new class of turbulence models is inherently 
grid- size-dependent, because increasing the grid resolution allows 
smaller and smaller turbulent scales to be resolved. 

We have introduced and implemented a novel two-stage 
procedure to efficiently estimate the level of scale resolution possible 
for a given flow on a given grid for partially averaged Navier-Stokes 
(PANS) and other hybrid models. This is a two-stage procedure. In 
the first stage, a RANS simulation with a standard turbulence model 
(STM) such as ks is used to produce an estimate of f k over the entire 
grid domain. In the second stage, we supply f k for the selective 
application of a hybrid turbulence model (HTM) such as the PANS 
formulation in regions in which the grid density is sufficient to 
resolve a portion or all of the large-scale flow structures. In the 
present implementation, f k is a function of length scale and grid size 
that represents a characteristic length- scale ratio. 

In the present paper, we selected the subsonic high- temperature jet 
flows to calibrate and validate the PANS approach. This 
implementation was a first step toward adding a variable-resolution 
turbulence model capability to CFD codes. The PAB3D code can 
now be used to refine the PANS formulation and to conduct 
validation computations using a variety of simple and complex flow 
physics problems. This approach needs to be calibrated, verified, and 
validated for a wide range of flow problems such as different 
temperature jet, cavity, and others. 
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